Quantum phase transition in the multimode Dicke model 
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An investigation of the quantum phase transition in both discrete and continuum field Dicke 
models is presented. A series of anticrossing features following the criticality is revealed in the band 
of the field modes. Critical exponents are calculated. We investigate the properties of a pairwise 
entanglement measured by a concurrence and obtain analytical results in the thermodynamic limit. 
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Introduction. An ensemble of two-level quantum sys- 
tems, each coupled to the common radiation field modes, 
the Dicke model, ^ was introduced initially to investigate 
superradiant emission — a dramatic increase in the rate of 
coherent spontaneous emission of an ensemble of atoms. 
The model has found a variety of applications in quantum 
optics, nanoscale solid-state physics, and quantum infor- 
mation theory. The critical behavior of this model was 
first discussed in Refs. where the single-mode Dicke 

model has been shown to admit a second-order classical 
phase transition. An investigation of the thermodynamic 
properties of the multimode model of superradiance was 
done in Refs. Handli where the conditions for the ther- 
modynamic instability were derived. In Refs. and 
Rziizewski et al. have demonstrated that the presence of 
the term in the interaction leads to the disappearance 
of the classical phase transition in the single- mode model. 
The multimode case has been addressed in Refs. H and 

m 

The quantum, or zero-temperature, phase transition 
reveals itself by the presence of nonanalyticity in the 
ground-state energy of a quantum system as some param- 
eter of the system — e.g., the strength of the interaction — 
is varied. In addition the energy gap — e.g., between 
the ground and first excited states — vanishes at the crit- 
ical point as a power-law function. This kind of phase 
transition is driven by quantum fluctuations and re- 
lated to qualitative changes in the ground state of the 
system. A single-mode Dicke model^ has been shown 
to be one of the examples of such behavior J^ii^iii Re- 
cently Emary and Brandesi^ have considered a single- 
mode Dicke model with no rotating-wave approximation 
and demonstrated that at the quantum phase transition 
the system changes from b eing qu asi-integrable to quan- 
tum chaotic. In Refs. Il6lll7l the universality prop- 
erties of the single-mode Dicke model have been ad- 
dressed; it was demonstrated that the finite-size scal- 
ing hypothesis^^ works well for the single-mode Dicke 
model — it pertains to the same universality class as the 
infinite-range anisotropic XY model in a transverse field. 
The divergent correlation length is typical for the quan- 



tum as well as classical transition at criticality. Quantum 
systems, however, acquire additional correlations that do 
not have a classical counterpart. These are entanglement. 
In connection with the quantum phase transition in the 
single-mode superradiance model, it has been studied in 
Refs.[liii and 113. 

In the present paper we study a general model^i^i^>i^ 
in which the atomic system is coupled to an arbitrary 
(discrete or continuous) set of radiation modes with ar- 
bitrary coupling constants. Many recent publications 
have emphasized the achievements in creating a nearly 
single-mode environment, such as, for instance, high- 
Q cavitiesi^ One should, however, note that the cavi- 
ties used currently have essentially a finite quality fac- 
tor. In fact it is rather a tough problem in many cir- 
cumstances to create an environment with a single field 
mode<^ Therefore the multimode generalization provides 
a more realistic scenario. In addition, significant progress 
in solid-state photonic crystals, systems of artificial 
atomS)^ and ion-trap structures'^ have allowed exper- 
imental investigation of a rich variety of atom-field cou- 
pling regimes, in which case the multimode Dicke model 
may be an exciting target. 

Multimode Dicke model. We consider a system of M 
atoms distributed over some finite volume. Each atom is 
regarded within the two-level approximation. The atoms 
are identical, having the same transition frequency wqi 
and coupled to the electromagnetic radiation field via 
a dipolc interaction with the coupling constants gj. for 
the jth atom and fcth mode. We also assume that the 
system resides in a volume small enough such that the 
electromagnetic field is nearly uniform within the atomic 
distribution. Therefore the atom-field coupling constants 
depend weakly on the atom's number — i.e., 5^ ~ gu. The 
overall Hamiltonian of the model is 
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We use units with h ^ 1. The atom-field coupling 
constants are proportional to the inverse square root 
of the atomic distribution volume, V ^ M. We ex- 
plicitly show this I/VM dependence in the interaction 
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term, whereas are finite and independent of M . Note 
that Eq.(IT]) has parity symmetry; i.e. it commutes with 

n^exp iTT (^j^ 4 flfc + Ej cr^) • 

The model given by Eq.([T]) incorporates chaotic dy- 
namics and is not integrable even for a single-mode field. 
However, it becomes approachable in the thermodynamic 
limit M — > GO for a large number of atoms. In this limit 
the solution of the single-mode model was carried out 
in Refs. and [Tsl where it was demonstrated that the 
system undergoes a quantum phase transition. Below 
we present a solution to the multimode model in a form 
convenient for further analysis of criticality and entan- 
glement. 

Effective Hamiltonian. Taking into account the struc- 
ture of the Hamiltonian ^ we introduce the collective 
spin operators^ J± = ^ ■ a^^ and Jz = \ Yli The cor- 
responding square of the total spin operator, J^, with 
eigenvalues j{j -f 1) commutes with and therefore 
j = Af/2 is a good quantum number. In our analysis 
we are interested in the situation when all atoms arc at 
the ground state before the transition, so that the ex- 
pectation of Jz , proportional to the magnetization of the 
atomic system, is —M/2. Using the Holstein-Primakoff 
transformation^'^ 



J+ = bWM-bfb, J-=JI, 
M 

T' 



J: 



(2) 



where b and 6^ are the usual bosonic operators, one can 
expand the Hamiltonian llj in powers of 1/M. Neglect- 
ing the terms of order 0(M^*), s > 1/2, we obtain the 
effective Hamiltonian in the form 

Heff=^uJob^b + ^ujkalak+^9k{ai + ak){b^ + b) -ujoj. 

k k 

(3) 

When the field modes are distributed continuously and 
separated from ujq by a gap, the first term represents a 
localized impurity. On the other hand, placing ujq within 
the field spectrum leads to the situation when confine- 
ment of the corresponding state may be defined only for 
a certain lifetime, represented by the imaginary part of 
the "energy." In the essentially discrete model of the 
field modes, k = 1,2, . . . , N, which the single-mode case 
is an extreme example of, making loq equal to one of the 
field modes energies ujk does not change the situation 
dramatically. Nevertheless, the diagonalization of Eq.Q 
can be carried out in the discrete case, provided proper 
integrations are introduced in the final equations when 
the continuous-model result is sought. 

The Hamiltonian ^ is quadratic and can be brought 
to the diagonal form H' — ^aqtla by the Bogolubov 
transformation defining a new set of bosonic variables 

13 a 

(4) 



where Aap, Aafs, Bpa, and Bpa are complex amplitudes 
representing the direct and inverse transformations. For 
convenience, we introduced Greek indexes a,(3,. . . , such 

that a^a=o = ^^^^ ^'^d a^^2:fc = ^fc^^ assuming that k, as a 
discrete index, counts the modes beginning with "1," as 
suggested earlier. 

The easiest way to proceed with the diagonalization 
is to use the commutation approach. One equates the 
coefficient before the creation (annihilation) operators in 
[aQ,77gyy] = f\ as soon as these commutators 

are evaluated using the direct (inverse) transformation 
(|4|) . This procedure leads to the system of equations 



+ '25akgkA+p + 25m^ 9k Alp = spA^^ , 

k 



+ 'A„n and A~ 



(5) 



where A'^p = Aap 



a/3 ana = A^p - A^fj. The 
system becomes complete provided the bosonic commu- 
tation relations \qa , 9^] = <5q/3 and , ajj] — Sa(3 are 
taken into account. One can also benefit significantly 
from the fact that B^fs = A^^ and Bafs = —Apa- This 
is easily obtained playing around with the commutations 
of Oq, and with Eq.Q. The complete solution of the 
system ([5]) is 
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The inverse transformation is obtained recalling that 
Baf3 = A*p^ and B^fs = -Ap^. 

Below the critical point. For gt > 0, the spectrum of 
the model is given by the solutions to the equation 
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For the single- mode case (iV =1) the solution coincides 
with one derived in Ref. Il3l . In general, for any N , the 
above equation cannot be solved analytically. However, 
an important property can be revealed from a mere anal- 
ysis of Eq. ([7]) . It can be noticed that the excitation ener- 
gies £ are real only for 7 < 7c, where 7 = J2k9k/^k ^^'^ 
7c = wo/4. This defines the critical point, with the re- 
lation between the coupling constants 4,^f.g'^/ujk = wq. 
This critical value of the interaction defines the quantum 
phase transition point, separating normal and superradi- 
ant phases. After the critical point 7 > 7c, the energies 
are found by properly displacing the bosonic modes, as 
will be shown shortly. For the moment let us discuss the 
continuum model corresponding to Eq.©. 

In the case of the continuum there is a difference be- 
tween the bosons created by the impurity operators b^ 
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and the field (bath) operators a|.. The bath energies are 
not affected significantly; i.e., one can reasonably assume 
that the impurity does not change the spectrum of the 
bath, LUk- At the same time, the bath does dress the 
impurity state. Further we assume that the impurity is 
localized and the corresponding dressed energy is outside 
of the field band. 

For the continuous field spectrum model, in Eq.([7]) one 
replaces ^ daj'D{uj) and gk g{'^)i where 

X>(a;) is the density of the bath modes. As an instructive 
example, we model the density of modes and the cou- 
pling constants as V^g^ (w) — auj'^9{uj — r2i)6'(il2 — i^) 
with n = 1, which corresponds to the well-known Ohmic 
case.^ The width of the band is, then, A = — f^i. 
This choice is for convenience only. It does not affect the 
qualitative structure of the result. Finally, the impurity 
level is given by the transcendental equation 

2aeu;o In + - + AauoA = . (8) 

Here the parameter a represents the strength of the cou- 
pling and can be varied. The critical point, as follows 
from Eq.®, is at = uJo/AA. Note that in Eq.® 
the bandwidth has to be finite for the quantum phase 
transition to appear — i.e., for ac > 0. In most phys- 
ical systems, however, the effective width of the band, 
Aefi, is finite due to the natural cutoff that comes from 
the density of modes, ^{uj), or the couphng constants 
g{cu) — i.e., duif{ui) ^ Aeff where f{uj) is some cutoff 
fimction. As a result, one simply has ac = loq/AAcs. 

Above the critical point. The quantum phase transition 
results in a qualitative change in the structure of corre- 
lations in the system's ground stateJ^ As the interaction 
strength g (7) or a increases beyond the critical point 
the parity symmetry of the system is broken and all the 
oscillators in Eqs. ([1]) and ^ obtain new equilibrium 
positions; i.e., the operators and a]j acquire c- number 
shifts. 

To get the effective Hamiltonian above the critical 
point we displace all the bosonic modes in Eqs. H]) and 
([2]) by Ofc ^ Cfe + \/jak and b ^ b + where Uk and 
(3 are some complex constants and the factor ^/j is for 
convenience of notation. One can easily show, following 
the logic of Ref. [13, that au and (3 are 0(1) quantities, 

/? = - 70/7, afc = - . (9) 

The effective Hamiltonian above the phase transition be- 
comes 

^eff^ ^^^ka\ak +ujQb'^b + ^^gk{a\ + ak){b + b'^) 

k k 

+ C(6' + 6^') + e, (10) 



where 

cjo (5 + 277 ri^) ~ 

4,(1+,) ' ^^•-'^Vtt^^'^' 

c^o (1 - 77) (3 + ry) 
877(1 + ,) ' 

^ 877(1 + 77) 27/ • ^ ' 

Here we also introduced 77 = "fc/^, recalling the earlier 
definitions 7 ~ '^kdtl^k and 7c = ^0/4. At the critical 
point, when 77 = 1 this effective Hamiltonian coincides 
with Eq.®. 

The diagonalization procedure here is similar to the 
one below the critical point. One makes the linear 
transformation of the form fl]) which leads to H'^^^ — 

^alliQa- The spcctrum above the critical point is de- 
fined by the solution to the equation 

2 ^ 2 

ii^V4^+e^-167' = 0. (12) 

At the critical point, when 7 = 7c(= luq/A), this equation 
coincides with Eq.©. The coefficients of the transfor- 
mation ([¥]) above the critical point can be obtained from 
Eq.® replacing: uq ^ ujq - 2C, gk gk- 

Equation ([7]) together with Eq. gives the complete 
spectrum of the multimode Dicke model in the quasi- 
classical limit for both normal and superradiant phases. 
Figure 1 illustrates the behavior of the energy levels be- 
fore and after the phase transition for a finite number of 
modes and gk = gy^ujk/2, so that 7 = Ng'^/2. The latter 
choice is for convenience only. One observes that the ex- 
citation energy originating from the atomic level hits the 
ground state at the critical point, which results in the 
divergence of the spatial parameters, as will be demon- 
strated shortly. After the critical point it goes trough the 
discrete "band" of the field-originated energy levels, cre- 
ating an anticrossing structure with each of them. Finally 
it formes a detached state developing the "gap" above the 
field levels. The energy gap for each anticrossing feature 
is inversely related to the number of field modes. The 
energies experience a jump in the derivatives at the crit- 
icality. 

It is interesting to note that for large gk, the creation 
and annihilation operators qn and (/jy corresponding to 
the separated top energy level become equal to a combi- 
nation of flo and aj only — namely, uq {3qN — q]^)/\/8. 
At the same time qa<N and ql^^j^ converge to —ak and 
— a|,, such that with the corresponding levels they resem- 
ble the initial {gk — 0) set of field excitations up to the 
phase factor. One concludes that at large gk the atomic 
subsystem and field states are factorized in each eigen- 
state of the system. The atom-field interaction merely 
dresses the atomic ensemble with energy s Aj. There- 
fore, one would expect zero entanglement between the 



4 



atoms and the field by that point. The anticrossing fea- 
tures following the phase transition can, then, be easily 
understood in terms of the transfer of an excitation while 
the coupling constants are changed adiabatically. For in- 
stance, if one has the top field mode excited initially, 
the adiabatic increase of the coupling will result in the 
transfer of this excitation to the atomic system with the 
energy above the field "band" and vice versa; see Fig. 1. 
Note, however, that all operators aa<N and aj^^^ are 
shifted after the critical point. 

Below the critical point the field is empty, (a[.afc) — 0, 
and the atomic system is in its ground state, {Jz)/j = 
— 1. Above the critical point the atomic system loses the 
magnetization as {Jz)/j = ~"fc/"f and the field becomes 
occupied as (a[.afc) — juso (7^ — 7^) /277c. This is similar 
to the single- mode case."'^'^ 

Let us now proceed to the continuous spectrum with 
the Ohmic bath model above criticality. Similarly to 
Eq.® we find 

The complete energy spectrum of the continuum-field 
case is presented in Fig. 2. It consists of two branches 
located below and above the energy band. The effec- 
tive impurity level vanishes at ac- As the strength of 
the interaction increases the dressed level approaches the 
continuum. At the same time a new dressed impurity 
state splits off from the field continuum developing the 
gap above the band, e 4aA. This can be understood 
comparing the discrete and continuum spectrum model; 
cf. Figs. 1 and 2. For larger bands — i.e., larger A — the 
split-off energy converges to the limiting linear behavior 
faster. 

Placing the impurity above the field band, one observes 
similar behavior; see Fig. 2. The over-band impurity level 
in this case diverges from the band starting with a = 0. 
The effective impurity level below the band splits off the 
field modes and vanishes at the critical point a — ac- 
Then it returns back to the field band. For both positions 
of the undressed impurity level the slope of the excitation 
energy above the band jumps at the criticality. 

When the field band is not bounded from above, except 
for the gradual effective cutoff, the impurity excitation 
level placed in the gap below the field band still undergoes 
the transition at the critical point, which is defined as 
ac — CLio/4Aeff, and then merges with the field band. 
The split-off level above the band is not developed in 
this case. 

Critical exponents. An important property of phase 
transitions is that systems with different microscopic dy- 
namics behave equivalently at criticality. Their behavior 
depends only on the dimension of the system and the 
symmetry of the order parameter j^i This phenomenon is 
known as universality. The parameters describing this 
critical behavior are universal quantities. Their behavior 
at the criticality is completely described in terms of a 
critical exponent. 
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FIG. 1; Energ y spectrum for the discrete model; A'^ = 10 and 
gk = g^/uJkJ^- The top-left inset shows the ground-state en- 
ergy of the system. The bottom-right inset presents the area 
around the critical point, in which the first excited state col- 
lapses to the ground state. It also magnifies the anticrossing 
between the effective impurity level and the first field mode. 
In all the frames the vertical dashed line shows the quantum 
phase transition point. 

Let us analyze the behavior of the lowest excited state 
energy (relatively to the ground state) in the vicinity of 
the critical point. We consider the discrete case first, fac- 
toring the magnitude of the interaction for convenience 
as gk = gXk- Here g is varied in the vicinity of the 
critical value gc = ^/ujq/ ^4^^, xi/^k- This is feasible 
as far as the shape of the dependence stays the same. 
Stepping away form the critical point infinitesimally as 
g — g^i: Sg, where the negative sign is taken for Eq.([7]) 
and the positive for Eg. p^ . one obtains 

e^Nclg-g^r. (14) 

Here zv — 1/2 is the critical exponent. The nonuniversal 
constant of proportionality which defines the energy scale 
from the left of the critical point is 

\gc + 4^oglEk=ixl^k J 

After the critical point it is greater by a factor of \/2. 
This behavior marks the second-order quantum phase 
transition. 

In the case of continuously distributed field modes with 
the Ohmic model one varies a around the critical value 
ac ^ ojq/ 4A in Eq. ([8]) and (fTS]) , observing similar behav- 
ior for the impurity excitation energy 

e^N^\a~acr, (16) 

with the same critical exponents z = 2 and ly — 1/4. The 
nonuniversal coefficient iV^ to the left of the critical point 
is 
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FIG. 2: Energy spectrum for the continuous field distribution. 
The field band is shown by the shaded area. The solid curve 
represents the effective impurity excitation energy. It begins 
below the band at a = and ends above the band at a ^ oo 
as purely atomic excitation energy. The quantum phase tran- 
sition occurs when the first excitation energy of the system 
(solid bottom curve) hits the ground state at Oc = A/4. The 
ground-state energy changes similarly to the discrete case and 
is not shown here. The parameters are ujq = A, ili — 2 A, 
and = 3A. The dashed top curve shows the effective im- 
purity when it is placed at luq = 4 A. In this case the quantum 
phase transition is manifested by the first exited level (bot- 
tom dashed curve) as it splits off below the band and hits the 
ground-state energy at Oc ~ A. 



After the critical point it acquires a factor of \/2. 

Entanglement as an order parameter. T\ie quantum 
phase transition is governed by a divergent correlation 
length at the critical point. This is similar to what hap- 
pens in the classical case. In the quantum case, though, 
there are correlations which do not have a classical ana- 
log. These purely quantum correlations are known as 
entanglement. It was demonstrate d^^i^^i^^ that there is 
a close relation between quantum phase transitions and 
entanglement. For a single- mode model it was shown^^ 
that the atoms become strongly entangled in the vicinity 
of the critical point. Osterloh et al.^^ have found that 
the concurrence is universal for the XY models. In what 
follows, we investigate the entanglement between the im- 
purity atoms in the multimodc field model and analyze 
how it is affected by the number of modes involved. 

The entanglement of an effectively bipartite system 
can be measured by the concurrence,^^ which is a convex 
function of the entanglement of formation. For a given 
two-atom density matrix it is defined'^^ as 



Cip) 



; |o, 2 max Ai 



(18) 



where Xj are the eigenvalues of the Hermitian matrix 
R = \/ \fPPsfP Etnd p = (ay®ay)p*{ay®ay). Due 
to the symmetry of the atomic system with respect to 
the exchange of atoms, one can evaluate pairwise con- 
currence, extracting two-particle states from the mul- 
tiparticle symmetric states. Wang and M0lmer have 



demonstrated-1 that the two-atom reduced density ma- 
trix written in the standard basis { |tT) , ITi) , lit) , lii) } 
takes the form 



u+ X 



+ 



Pii 



x\ u 

x+ w y* x*_ 

x+ y w x*_ 

U X- X- V- 



(19) 



As mentioned earlier, utilizing the symmetry of the 
atomic state under the exchange of particles one can ex- 
press all the elements in Eg. p^ — i.e., the expectation 
values {(Ju(J2j), etc. — via the expectation values of the 
collective spin operators;^ 



x±-- 



AP - 2M + 4 ( J2) ± 4 ( J^) {M - 1) 
4Af (Af - 1) '■ 
(Af-l)(J+)±([J+,J,]^) 



(A) 



A/2 



2M (Af - 1) 



4Af (Af - 1) 



y 



Af (Af 
J2) - Af 



2M{M - 1) 



(20) 



To calculate the above expectation values below and 
above the phase transition, we utilize Eq.Q in the limit 
of large Af . Below the critical point we obtain 



(Jl) = (Af - Vm) (b^b^) + 0{M-') 



(b^bb^b) + 0{M-''), 
+ 0{M 

+ {b^b) +0{M-'), s>l/2. (21) 



The other expectation values are either trivial (= 0) or 
can be derived from Eq. ((2T|) . After the phase transition 
we have 



(Jf)' = ^^^ + f (l-^) + M(l-ry)(6t2)' 

+ Af (1 - 277) (b^b)' + (b^'bh^'b)' + OiM-"), 



Af 



r;2) + M(,,_l) + ^(3^_l)(5t2)' 



+ ^(r/-l)(&^6)' + 0(Af-^), 
W = - + (b^b)' + 0{M-% s > 1/2. 



(22) 



We point out that the ground states in Eqs. (PT|) and 
arc different, which is the reason for the "primes" in 
Eas. (|22)) . Note also that after the critical point one has 
to include the shift of the bosonic operators in Eqs.Q to 
obtain Eqs. ([22]) . 

The expectation values of the bosonic operators be- 
fore the critical point are found utilizing the transfor- 
mation (HI) as (fotfo) = E/3^0/3> (&^^) = E/3 ^0/3^0/3, 

and (b^bb^b) = ^0/3^0/3^07^07 + E/3,7 ^0/3^07 + 

2 2 

E/3 7 ^0/3^07- Here we omit the complex conjugation 
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FIG. 3: Concurrence for the discrete model with different 
numbers of field modes. From right to left, A'' = 1, 2, . . . , 10. 
The other parameters are the same as in Fig. f . 



since all the amplitudes are real. Together with Eqs. © 
and ([7]) these relations give the expectation values of the 
total spin operators. A similar procedure leads to the 
explicit expressions for the expectation values after the 
critical point (P^ . 

For the Dicke model the reduced density matrix (fTO|) 
simplifies, x+ = x*^ = .t_ = x*_ ^ 0, y = y* — w, and 
the eigenvalues of the matrix R are Ai = 2w, A2,3 
± ^v+V- 1 , and A4 — 0. One can verify that A2 > 
Ai, w < ^/vTvZ, and 



\u ± ^v+V-\, and A4 — 0. 

A3 > A4 > Ai, w < y/vyvZ, and 2w < u + ,Jvyvl 
As a result the concurrence below and above the phase 
transition has the form C = 2(u — w). However, in the 
thermodynamic limit M ^ 00, the pairwise concurrence 
vanishes as C = 0{M~^) due to exchange symmetry, so 
we define the scaled concurrencei2i22iiH 



CMip) = 2M{u~w). 



(23) 



Using relations ((20)) . pTjl . and (|22p and the expressions 
for (b^b), (b^^), and (b'^bb'^b) we obtain the limiting value 
of the scaled concurrence. Coo = liniM-*oo C!m{p), hi the 
form 



where 



i^ + e) 



/3 



0- 



1 



(24) 



Wb = 



1+E 



(25) 



Here 9 = 1 below the critical point, 7 < 7c, and 9 = rj 
above the critical point, 7 > 7c. Also we imply the re- 
placements ujo ^ — 2C, gk — > gk for 7 > 7c. In 
Fig. 3 we plot Coo as a function of the coupling strength 
for various N, using the coupling constants in the form 
gk — gy^uJk/2 for an illustrative example. The concur- 
rence reaches the maximum and breaks at the critical 
point. We note that the increase in the number of the 



field modes leads to the corresponding increase of the 
maximum bipartite entanglement in the atomic system. 

Let us investigate the critical behavior of the concur- 
rence. Introducing an infinitesimal deviation of the cou- 
pling from the critical point g = gc i 6g one obtains 
(5 Coo = —WqSe/ujo from the first term in Eq. ([M|) . Here 
Wq is Wo, [see Eg. ipS]) ]. evaluated at the critical point. 
As a result one has 



I 



Coo W^N, 



9c 



(26) 



with the critical exponent zv — 1/2, which defines entan- 
glement (via concurrence) as an order parameter. One 
can easily generalize these results to the case of the con- 
tinuum field model. 

A^-term. In many circumstances the overall Hamil- 
tonian of the ensemble of two-state systems interacting 
with electromagnetic field would involve an additional 
field termr2i One would expect a Hamiltonian of the form 

H{t) = E-ii[P. - f A(t)]V2™ + E*ii where p, 

and Vj are momentum and coordinate operators of the 
jth subsystem with confinement potential Vj. Rewrit- 
ing the Hamiltonian in the eigenspace of the two-state 
systems with the assumptions discussed above we arrive 



— at Eq.(IT]) with an additional term M 



In terms 



of the standard bosonic operators it is Sk{a\. + '^fc)^ 
where Sk depends on the density of states and therefore 
is finite when M — > cx). Often this term is small with 
respect to the linear one and therefore is not consid- 
ered. Nevertheless, it was shownii^i^iiS that it can lead to 
the absence of the classical phase transition. Below we 
demonstrate that the quantum phase transition in the 
Dicke model with the term is present and can still be 
qualitatively described by the above results. 

We consider a single-mode model with loq — oJk = i^- 
The diagonalization of Eq.([T|) with the term — i.e., 
s{a^ + a)^ — yields an expression for the spectrum in the 
2s± v/(T 



form £2 = 1 + 2s ± v/(l + 2s^2 
e, s, and g are in units of uj 



- (1 -I- 4s - 4.g2). Here 
Exactly as before, one 



finds the critical point at g^ = VI + 4s/2, where the 
energy gap between the first excited and the ground 
state energy levels vanishes. Above this point the sym- 
metry breaks. At the criticality the energy scales as 
e — > \/45c/(l + 2s)|g — 3c|^^^- After the critical point 
there is additional factor of \/2 on the right-hand side. As 
a result, one observes the same transition with the same 
critical exponents and corrected position of the critical 
point. 

In the multimode case the derivation is more cumber- 
some due to the cross terms in . Considering the above 
strategy it is expected that the generalization will yield 
similar unimportant corrections. We leave a thorough 
proof of this statement to be given elsewhere. 

In summary, we investigated the model of the atomic 
ensemble interacting with an arbitrary discrete or con- 
tinuous set of bosonic modes. The system was shown to 
undergo a quantum phase transition of the second or- 



7 



der. In addition a series of anticrossing features follow- 
ing the criticality was revealed in the field band. Crit- 
ical behavior of the parameters was found explicitly. It 
was demonstrated that the pairwise concurrence of the 
atomic system is broken at criticality with the same crit- 
ical exponent as the energy. 



hospitability of the Center for Nonlinear Studies at Los 
Alamos National Laboratory. The work was supported 
in part by the NSF under Grant No. DMR-0121146. 



Acknowledgments 

The authors acknowledge contributing discussions with 
V. Privman and D. Mozyrsky. D.T. acknowledges the 



1 R. H. Dickc, Phys. Rev. 93, 99 (1954). 

2 K. Hcpp and E. H. Licb, Ann. Phys. (N.Y.) 76, 360 (1973). 
^ F. T. Hioc, Phys. Rev. A 8, 1440 (1973). 

* G. Comer Duncan, Phys. Rev. A 9, 418 (1974). 

^ K. Hepp and E. H. Lieb, Phys. Rev. A 8, 2517 (1973). 
® V. I. Emcljanov and Yu. L. Klimontovich, Phys. Lett. 59, 
366 (1976). 

K. Rztizcwski, K. Wodkiewicz, and W. Zakowicz, Phys. 
Rev. Lett. 35, 1967 (1976). 

* K. Rzgzewski and K. Wodkiewicz, Phys. Rev. A 13, 432 
(1975). 

^ B. V. Thompson, J. Phys. A 10, 89 (1977). 

^° B. V. Thompson, J. Phys. A 10, L179 (1977). 

S. Sachdcv, Quantum Phase Transitions (Cambridge Uni- 
versity Press, Cambrdgc, England, 1999). 

12 Y. K. Wang and F. T. Hioe, Phys. Rev. A 7, 831 (1973). 

1^ C. Emary and T. Brandes, Phys. Rev. Lett. 90, 044101 
(2003). 

M. Hillcry and L. D. Mlodinow, Phys. Rev. A 31, 797 
(1985). 

1^ J. Reslen, L. Quiroga, and N. F. Johnson, Europhys. Lett. 
69, 8 (2005). 

1® J. Vidal and S. Dusuel, Europhys. Lett. 74, 817 (2006). 
^'^ G. Liberti, F. Plastina and F. Piperno, Phys. Rev. A 74, 
022324 (2006). 

V. Privman, in Finite Size Scaling and Numerical Simula- 
tion of Statistical Systems, edited by V. Privman (World 
Scientific, Singapore, 1990), p. 1. 

N. Lambert, C. Emary, and T. Brandes, Phys. Rev. Lett. 
92, 073602 (2004). 



2° H. W. Chan, A. T. Black, and V. Vuletic, Phys. Rev. Lett. 
90, 063003 (2003). 

21 S.-Y. Lin, J. G. Fleming, M. M. Sigalas, R. Biswas, and 
K. M. Ho, Phys. Rev. B 59, R15579 (1999). 

22 A. Kiraz, P. Michler, C. Becher, B. Gayral, A. Imamoglu, 
L. Zhang, E. Hu, W. V. Schoenfeld, P. M. PetrofF, Appl. 
Phys. Lett. 78, 3932 (2001). 

2^ M. Keller, B. Lange, K. Hayasaka, W. Lange, and H. 

Walther, Nature (London) 431, 1075 (2004). 
2* T. Holstein and H. Primakoff, Phys. Rev. 58, 1098 (1949). 
25 A. J. Leggett, S. Chakravarty, A. T. Dorsey, M. P. A. 

Fisher, A. Garg, and W. Zwerger, Rev. Mod. Phys. 59, 

1 (1987). 

A. Osterloh, L. Amico, G. Falci, and R. Fazio, Nature 
(London) 416, 608 (2002). 
2^ T. J. Osborne and M. A. Nielsen, Phys. Rev. A 66, 032110 
(2002). 

2* L.-A. Wu, M. S. Sarandy, and D. A. Lidar, Phys. Rev. 

Lett. 93, 250404 (2004). 
2^ W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998); S. Hill 

and W. K. Wootters, ihid. 78, 5022 (1997). 
^° C. H. Bennett, D. P. DiViiicenzo, J. A. Smolin, and W. K. 

Wootters, Phys. Rev. A 54, 3824 (1996). 

31 X. Wang and K. M0lmer, Eur. Phys. J. D 18, 385 (2002). 

32 J. Vidal, G. Palacios, and R. Mosseri, Phys. Rev. A 69, 
022107 (2004). 

33 J. Vidal, Phys. Rev. A 73, 062318 (2006). 

3* The corrections for phonon fields may be more compli- 
cated. 



